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Abstract. Building upon previous results in scalar field theory, a formalism 
is developed that uses generalized Killing fields to understand the behavior of 
extended charges interacting with their own electromagnetic fields. New notions of 
effective linear and angular momenta are identified, and their evolution equations 
are derived exactly in arbitrary (but fixed) curved spacetimes. A slightly modified 
form of the Detweiler- Whiting axiom that a charge's motion should only be 
influenced by the so-called "regular" component of its self-field is shown to follow 
very easily. It is exact in some interesting cases, and approximate in most others. 
Explicit equations describing the center-of-mass motion, spin angular momentum, 
and changes in mass of a small charge are also derived in a particular limit. The 
chosen approximations - although standard - incorporate dipole and spin forces 
that do not appear in the traditional Abraham-Lorentz-Dirac or Dewitt-Brehme 
equations. They have, however, been previously identified in the test body limit. 



PACS numbers: 04.25.-g, 04.40.Nr 



1. Introduction 

Interactions between compact classical charge distributions and their own 
electromagnetic fields have now been studied in various contexts for more than a 
century [fl[l[lll[51[i[71[Hl[a[ini[ni[Il[Il[Il[IS]- The underlying physics is 
well-understood. Maxwell's equations are to be coupled to some version of continuum 
mechanics describing the type of matter under consideration. The resulting partial 
differential equations are then used to evolve the system forward in time from some 
appropriate set of initial data. Unfortunately, this proposal almost always requires the 
use of numerical methods. These are both difficult to implement and obscure general 
trends that may exist across classes of solutions. 

A trivial example of such a trend is found in the motion of sufficiently small 
charge distributions in the test body limit. To an excellent degree of precision, the 
center-of-mass position of such a body is governed by a single ordinary differential 
equation depending on only one parameter - the (conserved) charge-to-mass ratio. 
No other details of the object's composition need to be taken into account in order to 
characterize its bulk motion. This is not to say that they are completely irrelevant. 
Even in this highly restricted regime, the internal dynamics are not in any sense 
universal. Regardless, focusing attention on special variables such as the center-of- 
mass position can lead to interesting and relatively simple analytic results. 
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The example just discussed may be extended in several ways. Consider relaxing 
the constraints on a body's physical size (with respect to the wavelengths of externally- 
imposed fields) while remaining in the test body regime. Doing so, one finds 
excellent approximations for the motion - coupled to equations also describing the 
evolution of the mass and spin angular momentum - as ordinary differential equations 
involving a small number of parameters. These are the multipole moments of the 
charge distribution. The full theory of such effects has been developed by Dixon 
[TBI Ull US 033 ■ K holds even in strongly curved spacetimes, and is adequate for 
describing almost all isolated charges of practical interest. 

Still, there are some cases where a body's local field must be taken into account. 
Self-force or radiation reaction effects usually become significant only when changes 
in a charge's structure or overall motion occur on timescales approaching its light- 
crossing time. This might be expected in the highly-magnetized regions surrounding 
neutron stars and other extreme phenomena. Also of astrophysical interest is the 
gravitational analog of this problem. Consider, for example, a stellar-sized compact 
mass spiralling into a supermassive black hole. The gravitational radiation emitted 
in this process is expected to be observationally significant in the near future, so an 
improved theoretical understanding is vital. Considerable success has already been 
achieved in this area [101 [2Ql [21] , although some aspects might be better understood 
by first considering improved electromagnetic models. 

Regardless of motivation, this paper considers the problem of a compact charge 
distribution interacting with its own electromagnetic field (as well as those of distant 
sources). Despite the extensive literature on the subject, there are lingering questions. 
One of these regards the surprising success of point particle methods. It is commonly 
known that there are no solutions to the equations coupling Maxwell's equations 
to sources with vanishing radius but finite mass and charge. Despite this, various 
methods have been proposed to remove problematic infinities and obtain equations of 
motion that agree with those derived by considering well-defined limits of extended 
charge distributions [3 [9] . While there are intuitive arguments supporting these 
prescriptions, their connection to the full theory (involving extended charges) is 
unclear. The precise assumptions required to recover a physically and mathematically 
sensible "point particle limit" are not obvious, and should be better understood. 

The fact that point particle methods work at all is suggestive, however. Direct 
treatments of the dynamics of extended bodies have usually involved extremely tedious 
calculations in perturbation theory [21 [TTJ [12] . This has been especially true when 
allowing for arbitrary shapes and internal dynamics. The calculations used to obtain 
equations of motion for point particles are considerably simpler. This suggests that at 
least some of their hypotheses might have analogs in the theory of extended charges. 

We demonstrate this for the so-called Detweiler- Whiting axiom. Drawing on ideas 
originally developed by Dirac in flat spacetime [3], it states that a point charge only 
moves in response to a particular component of its self-field satisfying the vacuum 
Maxwell equations [9 . Applying this hypothesis very quickly leads to the Dewitt- 
Brehme (or generalized Abraham-Lorentz-Dirac) equation [4] typically expected to 
describe the motion of small charges in a fixed background spacetime. We show that 
a simple non-perturbative calculation allows an analogous statement to be derived 
from the full continuum theory. Furthermore, the "ignored" portion of a charge's self- 
field is shown to act mainly as an effective shift in its linear and angular momenta. 

Although self-force problems are typically thought of as intrinsically perturbative, 
several results derived here are exact up to effects related to the charge's gravitational 
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self-field. They are not immediately sufficient to compute the center-of-mass behavior, 
however. The main difficulty lies in obtaining suitable explicit solutions to Maxwell's 
equations. Effects of the spacetime curvature can also be difficult to take into account 
for large bodies. An approximation procedure is therefore applied in order to obtain 
explicit equations of motion, ft shrinks the charge's size while maintaining its charge- 
to-mass ratio and relative self-energy. The shape and internal structure are otherwise 
arbitrary except for the restriction that internal timescales do not shorten as the 
radius is decreased. Under these assumptions, equations describing the evolution of 
the body's mass, angular momentum, and center-of-mass position are derived. They 
effectively add Dewitt-Brehme-type self-forces to effects already expected for a test 
body up to dipole order. 

The methods used here generalize those recently developed for understanding 
extended bodies coupled to Klein-Gordon fields 22J. Roughly speaking, linear and 
angular momenta are first defined in curved spacetime by introducing a f O-parameter 
family of vector fields that generalize the Killing symmetries of Minkowski spacetime 
[23] . These ideas are reviewed in section [2j Various contributions to the evolution of 
the body's momenta are then analyzed in section[3] A particular component of the self- 
field is shown to act almost entirely as an effective inertia. Its remaining effect involves 
changes in an underlying Green function under flows defined by the generalized Killing 
fields. This has an intuitive interpretation as a violation of Newton's third law. It is 
shown that there are many interesting cases where such violations remain negligible 
or even vanish entirely. Finally, section [5] derives equations of motion for a small 
slowly-varying charge in curved spacetime. 



2. Preliminaries 



If the only long-range dynamical field in a system is electromagnetic, matter 
distributions that couple to it might be modeled by an action 

S = S M + I (J a A a - ^-F ab F ab )dV. (1) 

J 167T 

The field strength F ab = 2V[ a A^ is to be derived from the potential A a in the usual 
way, while the matter action Sm is assumed to be independent of these variables. 
Subtle issues regarding the proper formulation of electrodynamics in permeable and 
polarizable media will not be considered in any detail here [MJ [551 ■ It will 
be assumed that all such effects can be modeled by decomposing the current J a into 
appropriate "microscopic" and "macroscopic" components. Taking {1} as a given, the 
standard field equation 

V b F ab = AnJ a (2) 

follows immediately after varying the action with respect to the potential. 

Diffcomorphism invariance further implies that the system's total stress-energy 
tensor is conserved. It is useful to define a stress-energy tensor T ab only for the charged 
body under consideration. Let 

rp _ ^ SS „EM /q\ 

lab — . T— TT -L a b i W 



where 



TT^^(Fa C F bc -\ gab F'F,, n (4) 
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is the field's stress-energy tensor. 5S/5g ab represents a functional derivative of the 
action with respect to the inverse metric g ab . The support of T a f, is assumed to 
be a timelike worldtube W that may be foliated with a set of compact spacelike 
hypersurfaces. Total stress-energy conservation now requires that the body's energy, 
momenta, and stress satisfy 



The right-hand side of this equation may be interpreted as the generic force density 
exerted by an electromagnetic field on matter with current density J a . 

In general, much more can be learned from Sm implicitly contains a number of 
dynamical variables describing the detailed state of the material under consideration. 
Their behaviors can be extracted by applying the least action principle to each in 
turn. It is then found that systems described by the given action obey a coupled set of 
(usually) nonlinear partial differential equations. There are very few interesting cases 
where a direct analytical solution is feasible to obtain. Furthermore, these equations 
and their solutions both depend strongly on the type of material under consideration 
- i.e., the detailed form of Sm- Put another way, the full equations of motion always 
imply stress-energy conservation. The converse is almost never true. Despite this, it 
is known that there are cases where the details of a material's internal structure do 
not strongly affect its motion as a whole. Certain degrees of freedom can be almost 
entirely constrained by universal conservation laws like (|5|). 

Such laws are related to the underlying geometry (and the internal gauge 
symmetry, which only implies charge conservation here). As is well-known, the 
continuous isometries of Minkowski spacetime imply natural definitions for conserved 
quantities typically referred to as the body's linear and angular momenta. Somewhat 
unconventionally, we choose to view these objects as maps from the space of Killing 
fields into R. This is useful in that both the linear and angular momenta may be 
considered simultaneously. All quantities of interest are also scalars rather than first 
or second rank tensor fields. 

This scheme may be extended to curved spacctimcs that do not admit any Killing 
vectors. Given some reference worldline T and an associated foliation S(s), it is always 
possible to introduce something that may be interpreted as a generalized Poincare 
group [23j . All of its elements may be represented by linear combinations of 10 
independent vector fields in the four spacetime dimensions assumed here. By analogy 
to the situation in Minkowski (or de Sitter) spacetime, denote a body's linear and 
angular momenta by a map Pg(s) : GP xl^l. Let 



where £ a S GP is a generalized Killing field (GKF) [T7HT3I53]. Other reasonable 
definitions for exist. Terms may be added to the right-hand side of (0 involving 
J a and F a b, for example [T71 [18]. This possibility is discussed in the appendix. 
Furthermore, an additional effective momentum due to the body's self-field will be 
derived below. should therefore be thought of more as an initial guess than a 
fundamental definition. 

% This is only possible in a neighborhood of I\ It requires that any point in W n be connected 
to r n B(s) by a single (spacelike) geodesic. This essentially means that each spatial section of the 
body's worldtube is assumed to be within a normal neighborhood of the "origin." Spacelike geodesies 
are not easily focused, so this is a very minor restriction. 



X7 b T ab = F ab J b . 



(5) 




(6) 
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Relatively few properties of the GKFs are needed here. We note only that any 
element of GP is uniquely fixed by specifying £° and V Q £{, = V[ a £ b ] at any point 
7 £ r on the reference worldline. The full generalized Killing fields depend linearly 
(and non-degenerately) on these data. They also satisfy 

Ct9ab\r = V a £ 4 g &c | r = 0. (7) 
Any exact Killing fields that may exist are necessarily GKFs. A precise definition for 
these objects and an extensive discussion of their properties may be found in [23) . 

The map is easily related to more standard definitions of linear and angular 
momenta. Let these objects be tensor fields p a (s) and S ab = S^(s) on T. Then 
suppose that 

vds) = (Pae + \s ab v [a t; b] ) l{s) , (8) 

where the point 7(5) = T n S(s). The momenta are uniquely defined by this equation 
if is known for all possible GKFs. They are exactly the expected tensor fields in 
flat spacetime. Similar equations are also useful in Newtonian physics [22] . 

3. Self-interaction 

The motion of a compact charge is influenced by external electromagnetic fields as well 
as its own. While it is very difficult to compute every detail of the resulting dynamics, 
the momenta are relatively simple to evolve. Consider changes in V^(s) between two 
times s = s± and s — S2 > s\. Using ([5]) and (JBJ), 

MS = - = / (CJ b F ab + l -T ab C i9ab )dV. (9) 

JQ(si,s 2 ) z 

Here, f2(s!,s 2 ) is defined to be that portion of the worldtube bounded by the 
hypersurfaces S(si) and S(s2). It is trivial to derive from this an expression for 
the instantaneous momentum change — dV^/ds, although the averaged form given 
here will be more useful for now. 

The term on the right-hand side of ^ involving C^g ab represents the covariant 
gravitational force and torque on a compact mass. It exists even in the absence of 
any electromagnetic fields. Assuming that the body does not significantly influence 
the metric, g ab can usually be assumed to vary slowly inside the worldtube in (say) 
an appropriate normal coordinate system. It is then possible to perform a Taylor 
expansion and evaluate the relevant integrals term by term. Given ([7]), it is clear 
that the lowest-order approximation for the gravitational momentum shift involves 
the quadrupole moment of T ab . Full multipole expansions have been given in [T51ll9j . 
Specific examples of the dynamics arising from taking into account only the dominant 
quadrupole term have recently been discussed in several contexts [28l [2£l [30] . 

The focus here will be on the electromagnetic term appearing in This may 
first be simplified by splitting up the field into portions sourced by the charge itself 
and by the external universe: 

F ab =F:f + F s a f. (10) 

The self-field must satisfy 

V b Flf = AnJ a . (11) 

The current density here is assumed to vanish outside of W. All contributions from 
distant charges are therefore ignored in this equation. They are included in (J5J), 
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however. Suppose that the self-field is defined by solving (fTTj) with retarded boundary 
conditions. It follows that 

F*f = 2 / W [a Glf b ,J b 'dV\ (12) 
iff 

where G r °^,(x, x') is the retarded Green function for the vector potential in some (as- 
yet unspecified) gauge. These definitions imply that the external field is a solution to 
the homogeneous Maxwell equations at least in a neighborhood of W. It will often 
vary slowly throughout the body, so the forces and torques that it exerts may be easily 
approximated using series that involve successively higher multipole moments of J a . 

The influence of F^f 1 cannot be evaluated so directly. It will vary on scales at 
least as small as the body's diameter, so multipole expansions are not immediately 
helpful. Despite this, certain limits are expected to exist where the self-force becomes 
proportional to the square of the total charge q. This simple result - involving only 
the monopole moment of the charge distribution - suggests that there is a sense in 
which a simple multipole expansion is indeed applicable. The absence of self-forces 
and self-torques in Newtonian electrostatics (without dielectrics) leads one to expect 
that only a small portion of the field sourced by a relativistic body will directly affect 
its momenta. It is then reasonable to hypothesize that the remainder might vary 
so slowly as to exert forces and torques depending mainly on the body's first few 
multipole moments. This idea has been formalized by the Detweiler- Whiting axiom 
[9], which states that the motion of a point particle is only affected by a particular 
component of its self-field satisfying the homogeneous Maxwell equations. 

One of the main goals of this work is to write the laws of motion for extended 
bodies in a form where such results are essentially self-evident in an appropriate limit. 
The 'negligible' portion of the self-field is identified and (largely) removed at the outset, 
leaving only a nearly-homogeneous piece that produces the dominant self-force. This 
program first requires splitting the self-field into two components: 

F s a f=F! b +F* (13) 

By convention, the quantities on the right-hand side are respectively referred to as the 
singular and regular components of the physical self-field. These names refer to the 
respective behaviors of the two fields when associated with a point particle. Charge 
distributions considered here are assumed to be smooth, so all of their associated 
fields remain bounded. F^ b is intended to include most of the detailed structure of 
the self-field while having little effect on the body's overall motion. It will contribute 
an effective inertia, but does not otherwise influence the momenta as strongly as F^ b . 
One might expect this component of the self-field to be approximately Coulombian in 
form, and to be 'bound' to the particle in some sense. 

The regular component of the self-field is to be constructed so as to vary slowly 
inside the body under reasonable conditions. We therefore suppose that it satisfies 
the homogeneous Maxwell equations. Let 

V b F* = 0; V b F Q s b = 47rJ a . (14) 

Also demand that, for example, F^ b — 2V[ a A^ for an appropriate vector potential 

A^. Suppose that all such potentials are to be derived from Green functions in the 
simplest possible manner; e.g. 

A* = I G^J a 'dV. (15) 
Jw 
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Unless otherwise noted, we adopt the Lorenz gauge throughout this paper. The Green 
functions are then known to satisfy 

V 6 V fc G£ a ? - R a b G ( J = -4ng aa ,S(x,x'), (16) 

where g aa '(x, x') is any two-point tensor field that reduces to the metric &-S X ^ X . For 
concreteness, let it be the parallel propagator. Note that vector potentials obtained by 
substituting a Lorenz-gauge Green function into (fT5|) are only solutions to Maxwell's 
equations when V a A a = 0. Starting from the identity 

V a [g a a ,5{x,x')] = -V a ,5{x,x'), (17) 
it may be shown that [4] 

V[ '(V GH)=0. (18) 
The gauge condition is then seen to follow as long as V a J a = 0. 



3.1. Detweiler- Whiting Green functions 

A specific form for the singular self-field is already expected from the Detweiler- 
Whiting axiom. The appropriate Green functions are defined in Lorenz gauge, so they 
satisfy equations like (|16|) . Several auxiliary conditions effectively fix the boundary 
conditions, and therefore the solutions. First, the regular and singular Green functions 
are to be related to the retarded Green function exactly as might be expected from 

G%,{x,J) = G* a ,(x,x') + G s aa ,(x,x'). (19) 

The left-hand side vanishes when the source point x' is not in the causal past of 
the field point x. We restrict the S-type Green function from having support in the 
timelike past or future of any field point. This motivates its identification as the 
propagator for a bound field able to freely transfer momentum to or from the matter. 
Computing the singular field at a given point only requires knowledge of the charge 
distribution a short finite time into the past (and future). 

The last condition placed on the singular Green function requires it to be 
symmetric in its arguments: 

G s aa ,(x,x') = G s a , a (x,x'). (20) 

This is necessary in order for there to be a sense in which individual components of 
the singular self- force approximately obey Newton's third law (see section 15721 . The 
given conditions uniquely determine G^ a , and G^ a , , at least in some finite region. 
They are referred to as the Detweiler- Whiting Green functions. In flat spacetime, 
G L> = (Git, + Gfj)/2 and G*. = {G%, - Gfj)/2. 

More generally, these distributions can be at least partially written down with 
the help of Synge's world function a(x,x'). This two-point scalar represents one-half 
of the squared geodesic distance between its arguments. The support and symmetry 
requirements placed on the S-type Green function guarantee that it has the form 

Gl a ,=& aa ,5{o-)+V aa ,Q(a), (21) 

where A aa i — A a > a and V aa ' = V a > a are appropriate bitensors to be determined by 
substitution into (fTl)]) . Q(a) is a unit step function that vanishes when a < (i.e., 
when x and x' are timelike-separated) . The singular portion of (|2"Tj) is not difficult to 
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express entirely in terms of a and the parallel propagator g a a i . First introduce the 
(scalarized) van Vleck determinant 

A( X ,x>) = - det[ -^^ X ' )] . (22) 

It may then be shown that [TU] 

A Qa , = ^aa'A 1/2 , (23) 

where, again, g aa i is the parallel propagator. The remaining "tail" portion of the 
Green function cannot be evaluated so easily. A separate analysis must be carried out 
in each spacetime. V aa ' very rarely vanishes, although there are interesting cases - 
conformally-flat spacetimes, for example - where its existence is purely a gauge effect 
[31] . This contrasts sharply with scalar field theory, where physically- relevant tails 
almost always exist. 



3.2. Newton's third law 

It has been emphasized that the singular self-field is to be constructed so as to have a 
minimal effect on the charge's motion. Taking cues from the non-relativistic theory, 
one might expect to arrange this by ensuring that forces generated by F^ b satisfy 
an approximate form of Newton's third law. This procedure is considerably more 
ambiguous than in the scalar case described in [22], although it still provides useful 
motivation. 

Consider the force or torque exerted by one portion of the body on another and 
vice versa. This requires several clarifications. First, forces and torques are typically 
represented as vectors and 2-forms, respectively. The force on (say) one element of 
charge therefore cannot be immediately compared with the force on another element. 
This is resolved in exactly the same way that the linear and angular momenta were 
absorbed into a set of scalars via ([6]) and {8}. A particular component of force or 
torque can be picked out by choosing an appropriate GKF. Differentiating ([8]) and 
using (|Tj) shows that 

dV^/ds =Vt = F a £ a + ^N ab X7 [a d b] , (24) 

where the net force F a and torque N ab — N[ ab ] are defined by 

F a =p a - ^S bc j d R bcd a ; N ab = S ab - 2 P[a % } . (25) 

There is no loss of generality in considering for all possible GKFs instead of F a 
and N ab themselves. The generalized force acting on a particular component of the 
body may therefore be identified with its contribution to changes in the generalized 
momentum V^. 

Let T^{x^x') denote just such a generalized force density acting in a 'direction' 
£ a on charge in a small spacetime region dV due the singular field sourced by charge 
in dV. Given ([5]), it might be reasonable to define this by 

Fl(x,x')dV ^CJ b F! b [dV'} 7 (26) 

where -F^JdV^'] provides some notion of the field arising from charge in dV^'. It is 
perhaps most intuitive to define the field sourced by charge in a spacetime volume A 

by 

^ a s b [A]=2^V [a GS ]Q ,J b 'd^, (27) 
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so FfJdV] = 2V [a Gf ]b , J b 'dV. Note that the effective source in (gTJ) is usually not 

conserved when W A. The associated field strength is therefore not a solution 
to Maxwell's equations with any physically-relevant charge distribution. The field 
F ab [dV] apparently associated with charge only at a single spacetime point is actually 
generated by an effective current density 

J a = J a 6(x,x>) + ^-W a (W b G s w )J b '. (28) 

Despite these interpretational difficulties, the given constructions are still suggestive. 
Varying over all possible GKFs in (j2l)|) allows ordinary forces and torques to be 
recovered at least up to total divergences. Including additional terms that take into 
account the inherent ambiguities does not significantly alter the discussion. 

In terms of the force density the influence of the S-type electromagnetic self- 
field in ([9]) on the body's motion is easily shown to have the form 

f dVeJ b F! b = l [ dv( [ dV'[ff(x,x')+J$(x',x)] 
Jn £ J a, \Jw 

+ [ dV'[fl(x,x')-fl(x',x)]). (29) 
Jw\n / 

The arguments of £1 = 0(si, S2) have been dropped here for simplicity. Note that the 
first pair of terms on the right-hand side essentially measure the failure of Newton's 
third law. Using PO)) and (|26[) shows that 

Fl(x,x')+fl(x',x) = J b J b 'C^G s bb , - W a (J a eJ b 'G s w ) 

-V a ,(J a 'e'j b G s bb ,) 7 (30) 

where Lie derivatives of two-point tensors are defined to act individually on both 
arguments. If this result is averaged over finite volumes, the contributions from the 
two total divergences turn into boundary integrals that are easily dealt with. The first 
term more directly measures the degree to which action-reaction pairs are unbalanced. 
It is all that appears in the scalar field theory [55] . Importantly, the properties of the 
generalized Killing fields ensure that it remains negligible in many interesting regimes. 
It is also clear that it vanishes entirely if £ a is exactly Killing. This is a very useful 
consequence of adopting Lorenz gauge. Other choices can lead to Green functions that 
are not invariant under any Killing fields that may exist. 

3. 3. Laws of motion 

It is evident from this discussion that a significant simplification may occur if the 
only direct influence of the self-field can be shown to involve C^G^ a ,. This actually 
occurs quite naturally. Returning to |29|) . the left-hand side may be interpreted as 
the integrated force or torque acting on the body due to F^ b between times s\ and 
s 2 . It would therefore be reasonable to expect the integrand on the right-hand side 
to have essentially the same form on each time slice. More precisely, the integrand of 
the outer integral (over CI) in (f2"9")l should not explicitly depend on si or S2. It does 
do so, however. This is a consequence of having used the elementary identity 

f dV [ dV'A(x,x')= [ dV [ dV'A(x,x'), (31) 
Jn Jw Jn Jw\n 

which holds for any sufficiently well-behaved A(x, x') — — A(x',x). 
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The s-dependent terms can be interpreted as effective inertias. To see this, 
temporarily suppose that S2 — s± is sufficiently large that any element of S(s 2 ) may 
be connected to any element of £(si) by a timelike curve. Roughly speaking, this 
means that we're considering time differences larger than the body's diameter. The 
assumed support properties of G^ a , then imply that the second line of (|29|) becomes 
a double integral over points in W within a light-crossing time of S(si) and E(s2). It 
does not depend in any way on the behavior of the body in the distant past or future, 
nor in the "middle" of tt. Furthermore, the contributions from regions near £(si) 
and S(s2) can be separated into a difference of two integrals with identical integrands 
(they differ only in the region of integration). Combining ([9]), ([29]) . and (f30|). 



where 



dV 



T ab £ i g ab + Z a J b (F£ t + F»i,) 



dv'rr c^G s aa , 



, (32) 



%(s) = l(f dVJ a £tA s a \E-]- J^_dVJ a C 5 A s a [J: + ^ + Jj b Afj a dS a . (33) 

Bracket notation used here is defined by analogy to (|27|) . The four-dimensional 
volumes £~(s) and £ + (s) respectively denote the past and future portions of the 
body's worldtube with respect to £(s). 

Although this result was most easily motivated by considering momentum shifts 
over long times, it is actually valid for any choices of s% and S2- It is therefore possible 
to write (1321) in the instantaneous form 



d 



dS c t c 



T ab C^g ab + S, a J\K 



TiCXt 

ab 



dvrr c^G s aa , 



(34) 



t a represents the time evolution vector field for the foliation {£}. This is a general law 
of motion for compact bodies interacting with electromagnetic fields. Similar methods 
have been used to obtain a nearly identical result in scalar field theory 122]. 



3-4- Self-momentum 

The manner in which £^ appears in (|34[) suggests that it may interpreted as an effective 
momentum associated with the S-type self-field. Another (more obvious) possibility 
can be obtained from an expression like the one used to define V^. Replacing T ab in 
([6]) by the electromagnetic stress-energy tensor Q associated only with F^ b , let 

Ui = hL 9Cd{F ^ Fbd ~ \9a b g fh F? f F* h )edS b . (35) 

Unlike £{, this quantity depends on the details of the field far outside of the body's 
worldtube. It therefore requires knowledge of the charge distribution in the distant 
past and future. It also isn't clear that the integral converges if E is extended to 
spacelike infinity (even discounting the difficulty in defining £ a at large radii). 

is therefore unsuitable for use as the effective self-momentum of a dynamical 
charge distributior|§|. Despite this, it is known to give reasonable results at least in the 

§ It is still tempting to try, however. Consider an analog of l|34p with replaced by hi^. For 
definiteness, let each E be bounded, so all integrals are manifestly convergent. The force due to 
C^G S , then disappears. It is replaced by a flux integral through 9S and an additional body force 
requiring knowledge of C^g a t outside of W. Both of these terms are difficult to control. 
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stationary case. For simplicity, consider a system in flat spacetime where all variables 
remain invariant with respect to a timelike Killing field t a = d/dt. Also suppose that 
each £ is a hypersurface of constant t. It may then be shown that 



if the potential is in Lorenz gauge. Integrating ()33() by parts recovers this expression 
plus 



£j = under the given assumptions. The proposed linear and angular momenta of 
a body's stationary electromagnetic field are exactly as expected. 

Returning to the general non-stationary case, this correspondence together with 
(|34p provides strong support for the definition of an effective momentum 



from which to compute a body's motion under the application of external fields. 
Deriving (say) a center-of-mass worldline from rather than is closely related 
to the use of effective masses in more standard treatments of radiation reaction (see, 
e -g-, [HI H2 [32] ) • Fully taking into account all of the self-field's momentum - rather 
than only its inertia - is essential for understanding the dynamics of general charge 
distributions. 

One reason for this is that the field's energy distribution may have a different 
center than the matter's. The orbital component of angular momentum effectively 
forms the first moment of the self-energy distribution. It is especially clear in ([36]) 
that it measures the displacement of the "center-of-charge" away from the origin where 
the associated boost-type Killing fields vanish. In general, mass centers computed from 
■p£ and are therefore not the same. Simple behavior can only be expected to follow 
from carefully-chosen definitions. 

Furthermore, typical mass renormalization techniques that have appeared in the 
self-force literature effectively assume that the linear momentum associated with the 
self-field is parallel to that of the matter distribution. There are many interesting 
systems where this is not the case. As an example, a simple calculation using (|36p 
shows that the 3-momentum associated with F^ b can be finite even in the stationary 
case. Despite the apparent lack of motion, the mechanical momentum needn't vanish 
either. These two facts are closely related. Stress-energy conservation can be used 
to show that the mechanical and field momenta are actually negatives of each other 
in these cases. Choosing any translational Killing field £ a transverse to E, V$ = 0. 
The total 3-momentum of the body as suggested by this formalism therefore has the 
expected limiting behavior in time- independent systems. Nonzero values for by 
itself have been referred to as "hidden momenta" in the literature [32] [33] . Failing to 
properly take these effects into account when deriving laws of motion could lead to very 
peculiar predictions. This appears to have been a major reason for the complicated 
equations of motion derived in |12j . 

Although avoids these problems, its definition is rather peculiar. Momenta 
are usually obtained by integrating physical quantities over an appropriate spacelike 
(or occasionally null) hypersurface. The initial guess ([6]) for a body's mechanical 
momentum fits this model. £^ does not. Its definition involves an additional integral 
forwards and backwards in time off of the fiducial hypersurface. Computing the 




(36) 




(37) 



(38) 
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momentum at a given time s therefore requires knowledge of the charge distribution 
in a four-dimensional volume centered on j(s) and with a radius of order the body's 
proper diameter D. If the system of interest does not change significantly on these 
timescales, it is possible to evaluate the time integrals in (|33p explicitly. One is then 
left with the expected 3-dimensional integrals. 

Another brief aside regards the apparent (though fictitious) gauge dependence 
of It was noted that (|36|) holds only in Lorenz gauge. For an arbitrary vector 
potential generating F^ b , there would be additional terms involving V a A a i . Of course, 
the potentials we have been using are fixed uniquely. There is therefore no gauge 
freedom at all that our expressions could possibly depend on. 

Still, it is interesting to consider "Detweiler- Whiting type" Green functions 
satisfying all of the same properties as the usual ones except for the Lorenz gauge 
condition. First introduce a two-point tensor field Xa{x,x') that vanishes whenever 
its arguments are timelike-separated. A new Green function G^ a , can then be obtained 
from the original via 

Gl a > = G S aa , +V a >Xa + V aX a> • (39) 

This is symmetric, has the same support restrictions as G^ a ,, and generates the same 
field strength F*L. Attempting to use it in (|3"3"|) would yield an apparent momentum 

£e = £t + \l &s a j a f &v'j a 'c aa ,. (40) 

There are no changes if ^7[ a ^^Xb] = 0. It is only in this highly restricted sense that 
our equations can be considered to hold for different Green functions. It would be 
aesthetically pleasing to do away with the Lorenz condition, although it is not obvious 
how to do so. Everything is well-defined, so there arc no fundamental problems. 

4. Generalized Detweiler-Whiting axiom 

It is now possible to derive a generalization of the Detweiler-Whiting (DW) axiom 
introduced in [9]. As remarked above, this postulates that point charges move as 
though they were test bodies interacting with 

F l :r = + fS (4i) 

via the Lorentz force law. This field is a solution to the homogeneous Maxwell 
equations. It is almost always smooth at the particle's location, so the DW axiom 
provides a well-defined prescription for computing trajectories in given external fields. 
Maxwell electrodynamics is not compatible with the existence of point particles, so 
this hypothesis cannot be derived in a strict sense. 

Physically, however, point particles are usually introduced in order to approximate 
the motion of sufficiently small extended charge distributions. The DW axiom can 
therefore be viewed as a statement on the limiting behavior of a certain class of 
nonsingular objects. Interpreted in this way, it is amenable to derivation, and the 
formalism developed here does so very efficiently. Let the generalized electromagnetic 
Detweiler-Whiting axiom be the statement that the effective momentum defined 
by ([6]), ([33]) and ([38]) evolves like that of an extended test charge with the same T Q t, 
and J a , but moving under the influence of F^ m rather than Fj^. We say that this 
is true (with respect to £ a ) to order if 

^ = jf dSct < (±T a %g ab + ?fl*r) + H . (42) 
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It is interesting only when ej can be shown to be negligible compared to other terms 
here. 

Given (|34p . the critical question therefore becomes the relative importance of 

f dS c t c [ dV'J a J a 'C^G s aa ,. (43) 

JT, JW 

Firstly, it is clear that C^G^ a , = for any exact Killing field £ a that may exist. It 
follows that = in these cases. The generalized DW axiom with respect to any 
Killing field is exact. This means that there is a sense in which a particular component 
of a body's momentum is unaffected by F^ b . In maximally symmetric spacetimes - 
Minkowski and de Sitter - all components of linear and angular momentum react 
only to the regular self-field. This result holds even for charge distributions that may 
be large and experience rapid internal oscillations. The only approximation used for 
this conclusion is the assumption that the body does not significantly influence the 
spacetime geometrjfjj. 

It is more interesting to ask what happens when no Killing fields exist, =/= 
in these cases, although it can still be negligible. Consider a one-parameter family of 
matter fields that shrink down to a (fixed) central worldline T without distorting their 
shape. Introduce an orthonormal tetrad e\{x) in a neighborhood of these bodies, and 
let 

J a (A; s, r) = \- 2 e a A (s, r)J A (X, s, r/\) (44) 
T ab (X; s, r) = \- 2 e\{s, r)e h B {s, r)f AB {\, s, r / X) (45) 

for all A > 0. (s, r) are Fermi normal coordinates with respect to T. J A and T AB 
are both assumed to be smooth in all of their arguments. These equations do not 
change in any essential way under smooth A-independent coordinate transformations. 
Families of this type are essentially the same as those considered in [T5] . Somewhat 
more general ones were discussed in [22] . 

The methods required to estimate (|43j) are very similar to those used for its scalar 
analog in [22 . There, it was shown that at fixed r/A and r'/A, Lie derivatives of the 
various bitensors in ([21]) with respect to any GKF £ a fall to zero at least as fast as 

C^(x n )\ a=Q ~ 0(A 4 ), £ e lnA(x, 7 )U=o~0(A 2 ). (46) 
near T. Similar methods can easily be used to show that |(7J) also implies that 

£ 5 3V| CT =o~0(A 2 ). (47) 
The notation here means that, for example, 

lim £*a(a;,7)/A 3 = 0. (48) 

Together, these results imply that (|43p is 0(A 3 ). It follows that the error term in (|4"2")) 
scales like 

e 6 ~ 0(A 3 ) (49) 

as A — > 0. The Lorentz force due to the external field scales as A 1 , while the force due 
to F^ b scales like A 2 . For sufficiently small A, both of these terms dominate over (|43[) . 
This shows that the generalized Detweiler- Whiting axiom holds in an interesting sense 
with respect to all GKFs and all charge distributions described by (|4"4")) and ([43[) . 



| If it did, exact Killing fields could only exist in stationary and/or axisymmetric systems. 
Momentum evolution is not particularly interesting in these cases. 
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5. Motion of small bodies 

It is often useful to simplify an object's motion by considering only the behavior 
of a single worldline intended to be representative of its overall motion. Several 
reasonable possibilities for such a mass center have been proposed. Their differences 
are all expected to be bounded by distances S~f ~ (spin angular momentum) / (mass) 
[34] , Bodies obeying standard energy conditions likely satisfy D > £7 [35], so all 
potential center-of-mass worldlines should be contained in the convex hull of W. 
There is therefore a sense in which any of them is "representative" of the bulk motion. 
Despite this, a surprisingly large number of potential definitions are excluded simply 
by demanding that the motion of an uncharged body in flat spacetime be described 
by a unique geodesic (given standard initial conditions). We adopt a definition - 
originally due to Tulczyjew [36 ] - that does satisfy this criterion [17j . It also uniquely 
specifies both the central worldline T and the foliation S(s) in reasonable systems [37]. 

Briefly, we let E(s) be formed by the intersection of all geodesies that remain 
orthogonal to p a - extracted from using ([5]) - as they pass through 7(s). The 
boost momentum is also assumed to vanish: 

(p a S ab )\ r = Q. (50) 

This equation is a direct analog of familiar non-relativistic center-of-mass definitions. 
In general, the resulting p a will not be tangent to T. Deviations will usually be very 
small, although not completely negligible. It is possible to algebraically solve for the 
center-of-mass velocity 7° in terms of q, p a , S ao , the fields, and the higher multipolc 
moments of T ah and J a . This was accomplished exactly in [38] in the absence of any 
electromagnetic fields. A generalization to the present case is not difficult, although 
the result is rather unwieldy. We will only be using an approximation valid to 0(A) 
for the one-parameter family of charge distributions considered in section [4] 

It is first convenient to change definitions of the momentum so that all quantities 
are referred to the map V'^ introduced in the appendbjj]. Hats and primes will be 
omitted for brevity, however. Now let the worldline parameter s be normalized such 
that 

i a Pa = -m, (51) 

where the mass is defined by letting 

p a = mn a (52) 

and n a n a = — 1. These equations imply that s is not precisely a proper time unless 
7° = n a . With these conventions, the momentum- velocity relation is easily derived 
using the methods of [38 . Substituting (|A.7|) , it is approximately given by 

mj a = p a + {qS a c F^ - 2mQ c [a F$)n h /m + 0(A 3 ). (53) 

Although one might expect p a to be exactly 7717°, it was already remarked in section 
13.41 that additional terms must be present. A first approximation for these "hidden 
momenta" appears explicitly here. Such effects are not intrinsically electromagnetic. 
Similar relations exist even if J a = (although they then become nontrivial only at 
higher orders in A). 

% This is not strictly necessary. Multipole expansions for the force and torque associated with V'^ 
have already been derived 17, 18, 38 . Being able to use these results simplifies the present derivation. 
There are also some technical properties of these definitions that are conceptually attractive. Still, 
exactly the same analysis could be carried out using the original (renormalizcd) momentum Vt- 
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An immediately useful consequence of (|53[) may be derived by contracting it with 
7°. This shows that 7 a 7 a = — 1 + 0(A 2 ). The worldline parameter s is therefore a 
proper time up to terms of order A 2 . These corrections are too small to be relevant to 
any results derived in this section. The distinction between s and a true proper time 
may therefore be ignored. 

In any case, the full equations of motion can now be written down very easily. 
One of these is (|53]l . The other two are obtained by combining (|25|) . (|A.6j) and (|A.7|) 
to yield 

Pa = (qF*r - \s cd R abcd )i b ~ \Q bc V a F°? + 0(A 3 ), (54) 
S ab = 2( P[a j b] + Q c [a F^f) + 0(A 3 ). (55) 

All three equations must be solved simultaneously in order to determine the motion 
(as well as the mass and angular momentum). They are not complete, however. 
Appropriate prescriptions must also be given for Q ab and F^ b . In any case, the only 
effect of the self-field is to act as a small time-dependent shift in the effective external 
field. The dynamics here are therefore very similar to those of a test body in the 
dipole approximation. 

A better understanding of this system may be gained by introducing (generalized) 
Fermi derivatives along T [IT]. For any vector field k a , let 

5^- = k a +2n^ a n h] k b . (56) 
as 

The Leibniz rule together with Dpg ab /ds — may used be to generalize this definition 
for tensor fields of arbitrary rank. Intuitively, it measures changes in components 
measured with respect to a non-rotating zero-momentum frame. 

The Fermi derivative of the angular momentum takes on a very simple form. 
Before presenting it, note that the center-of-mass condition allows S ab to be replaced 
by a spin vector 

S a = -\e abcd n b S cd . (57) 

Inverting, S ab = £abcdn c S d . A simple substitution now shows that 

^ = -e ab cd n b Q\ c F$ + (A 3 ). (58) 

The right-hand side represents the ordinary torque expected to act on a test body 
with electric and magnetic dipole moments. Spin changes absorbed into the definition 
of the derivative operator are due to Thomas precession. They are necessary in order 
to preserve p a S a — along T. 

Another interesting consequence of the equations of motion is an evolution 
equation for the mass, m is not usually conserved, but rather varies according to 

±(m- \Q ab F^ = -l^lFtf + 0(A 3 ). (59) 

One sees that a charge's electromagnetic potential energy directly contributes to its 
inertia. The right-hand of this equation is related to the internal mechanical work 
performed by an external field. Phenomenologically, one might consider a dipole 
moment whose components rigidly rotate with an angular velocity tied to S a via some 
(similarly rotating) inertia tensor. It may then be shown that m is given by the sum 
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of a constant "internal energy," the dipole potential energy, and the body's rotational 
kinetic energy [T?] . 

Perhaps the most interesting quantity to calculate is the center-of-mass worldlinc 
r. It is difficult to get an intuitive feel for solutions to the coupled system we have 
derived. Consider instead the derivative of (|53|) with appropriate substitutions from 
(|54| and (|55p . n a may be eliminated by replacing it with 7° + 0(A). Making this and 
related substitutions results in 



The first three terms here are well-known. They are respectively the Lorentz (with 
self- field), Papapetrou and electromagnetic dipole forces. It is unfortunately rare - 
though not unheard of [2H HE] EH] ~ to emphasize that such forces represent rates-of- 
change of momentum rather than direct influences on 7717°. Taking into account the 
aforementioned hidden momenta yields a number of additional terms. These are all 
test body effects that have nothing to do with self-interaction. All such terms were 
implicit in Dixon's work [17j . although they were not displayed in this way. Various 
other approaches have found similar effects [MJ [3H [Ml SOI HI] ; although there is a 
surprising amount of disagreement in the literature on relativistic dipole forces. Very 
recently, however, (|60[) was also derived (with radiation reaction) in flat spacetime 
using a completely independent method |15j . Some analysis of its consequences may 
be found there. 

It should be emphasized, however, that (|60p should not be used to indiscriminately 
replace ([55)1 and (jM)) . There is no guarantee that its solutions are physically 
meaningful over long times. They may not even exist. The implicit assumption 
that (say) 7° has a well-defined limit as A — + could also break down. At least 
if is ignored, no such problems arise with the original system. The initial 
value problem has fundamentally changed through the use of substitutions such as 
qF°£ t n b — ► mj a + 0(A 2 ). Similar problems do arise in the original equations of 
motion if the usual form for the regular self-field is used. This introduces another 
7 term multiplied by a small parameter; cf. (|61[) below. One can see that many of 
the well-known problems of the traditional Abraham-Lorentz-Dirac equation actually 
arise even in a slightly careless discussion of test particles with significant dipole 
moments. A simple recipe for avoiding these problems is to use an order reduction 
procedure. Justifications for this point of view as well as more rigorous applications 
of perturbation theory may be found in [15j and [41] . 

Putting aside technical issues, we have not yet discussed the self- force in any 
detail. As can be seen in both (|54|) and ([60]) . this is given by a simple Lorentz-type 
force associated with F^ ) ( r y(s)). The regular component of the self- field is known to 
have a well-defined limit even for a point source with finite charge. Explicitly, it is 
given by [TO] 



The existence of such a limit suggests that the regular self-field of a reasonable 
extended charge distribution varies negligibly on each time slice Y,(s)DW. We therefore 
conjecture that (fB"Tj) describes the field of an extended charge distribution up to 0(A 2 ). 



-(m 7 „) = (qF*r ~ ^Ra bcd S cd )j b - -Q bc V a Fi 
as 2 2 

- 2 \Q c [a F^ b + H£(Q c [a F^ 7 b )j + 0(A 3 ). 



nCXt 

be 



Sabl" 



(60) 




(61) 
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Substituting it into (|60p while reducing order and coupling (|55p describes the motion 
of small self-interacting charge distributions in fixed background spacetimes. 

6. Conclusions 

We have developed a general formalism that greatly simplifies the self-force and self- 
torque problems for compact charge distributions coupled to Maxwell fields. The 
central insights regard the effects of the (Detweiler- Whiting) "singular" component of 
a charge's self-field on its motion. As shown by there are two basic mechanisms. 
First, F^ b contributes an effective inertia to its source. This component of the field 
makes the body appear as though it had additional components of linear and angular 
momentum. These depend only on a quasilocal knowledge of the charge distribution, 
and cannot generally be mimicked by simple mass shifts. 

This is the only direct effect of the singular self-field in maximally symmetric 
spacetimes. More generally, though, it also exerts a force depending on Lie derivatives 
of the Green function G^ a , with respect to the generalized Killing fields. This has a 
simple interpretation in terms of the degree to which Newton's third law is violated 
by forces associated with F^ b . All of these conclusions are exact (up to gravitational 
backreaction effects), and were derived without any reference to specific solutions 
of Maxwell's equations. They apply for charge distributions of almost arbitrary 
generality. Only very weak restrictions need to be assumed regarding size, shape, 
self-energy and internal composition. 

Specializing to charge distributions that are small and have slow internal 
dynamics allowed the explicit (though approximate) equations of motion (|53|) . (|54[) 
and (|58p to be derived. As expected, these involve the DeWitt-Brehme self-force. 
Electromagnetic and gravitational dipole forces are of similar magnitude in the 
approximations used, so all effects associated with dipole test bodies are also recovered. 

The dynamics are "universal" in the sense that they depend on the internal 
structure only through a handful of simple parameters. A similar statement is unlikely 
to hold if the equations of motion are expanded to higher orders in A. Instead, 
coefficients would probably arise that cannot be adequately approximated by any 
finite number of multipole moments. The formalism developed here likely provides a 
useful starting point for studying such effects. Still, it is not clear how useful such a 
derivation would be. We have already noted that a number of reasonable definitions 
for a body's momentum and (and therefore) mass center may proposed. These could 
all satisfy Detweiler- Whiting axioms in the sense of section [U It is likely, however, 
that higher-order self-force and self-torque effects would depend on such choices. It is 
therefore unclear in what sense they would be physically relevant. 

Additional work is also required in order to understand gravitational self-forces 
in the context of the formalism presented here. This problem is important in its 
own right, although ignoring it has also constituted a significant gap in the equations 
of motion presented in section [5j Dimcnsionally, gravitational and electromagnetic 
self-forces should enter at the same order of approximation for families of matter 
distributions described by (|44|) and (j45j) . There could then be nontrivial couplings 
between the two types of self-interaction. As an example, recent work on the motion 
of a charged black hole suggests that the force coupling to the background Ricci tensor 
in (|6ip is modified when gravitational backreaction is taken into account [14] . It would 
be interesting to obtain a more local understanding of this effect. 
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Appendix: Alternative definitions for the bare momentum 

There is considerable ambiguity in defining the momentum of a charged particle. Many 
possibilities look at least superficially reasonable while still satisfying a generalized 
Dctweiler- Whiting axiom like the one discussed in section!?] In principle, each element 
of ([6]) may be criticized. Temporarily leaving aside the overall form of that equation, 
definitions for both the body's stress-energy tensor T ab and the generalized Killing 
fields £ a can easily be argued. The only properties of GKFs used here were that they 
satisfy ([7]) and form a ten-dimensional group which includes any exact Killing fields 
that may exist. These conditions may be satisfied by vector fields different from those 
considered in [23) , although any alternatives would remain similar on sufficiently small 
scales. 

We have defined the charge's stress-energy tensor by (|3]). Again, other reasonable 
possibilities exist. There is a sense in which only the full stress-energy tensor is 
physically relevant. Some interesting perturbative statements can be made using this 
total [TS], although more generally, some split between material and electromagnetic 
components is essential. The only obvious constraint on the body's stress-energy 
tensor is that it must vanish outside of W. While the electromagnetic stress-energy 
tensor is unambiguous in vacuum, it could reasonably be given a different form inside 
matter. It is unclear how the various possibilities would affect the forms of the 
equations derived in this paper. 

Separate from these issues is the basic structure of ©• It is already strongly 
suggested by the introduction of an effective momentum V% in section T3. 41 that V% is 
incomplete in some sense. A detailed study of the possible multipole moments that 
can be associated with T ab suggests that V% should be replaced with [17l [18] 



^(«)=P c («)+jf ^ K- x SJ a ,o{x' , 1 )^(x')F a ' b ' (x')d^j J a (x)dS a , (A.l) 

It is implicit here that 7 = j(s). x' = x'(k) is an affinely-parameterized geodesic 
satisfying x'(0) = 7 and x'(l) = x. In terms of a potential A a , VI takes on the more 



intuitive form 

V*=Ve 



L Jo 



J a dS a . (A.2) 



There is strong evidence that these definitions also arise naturally from Lagrangian 
considerations [42]. Note that the {£ a A a )\ x term also appears in the effective 
momentum £e defined by (|33j) . 



Although the relation between V'^ and the continuum variables T ab , J a , g a b, and 
A a is quite complicated, it has a number of useful properties. There is a simple relation 
to conservation laws, for example. Temporarily suppose that there exists a vector field 
ip a satisfying C^g a b = C^Fab — 0. It is then possible to choose a potential A' a such 
that £^A' a — 0. It immediately follows from ([9]) that 

^[^ + «)W = 0. (A.3) 

where q is the body's total charge. The quantity in brackets is conserved for all 
extended charge distributions. Different choices for A' a can shift it only by an overall 
constant [T7]- A special case of this result applicable for pointlike test particles is 
commonly derived in textbooks [32] , 
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Another advantage of V'^ is that its changes have a particularly simple form 
in terms of the multipole moments of T ab and J a [18] . Unfortunately, equivalent 
expressions are quite complicated when expressed directly in terms of these fields. 
There is no fundamental obstacle to using (|A.1[) instead of ^ in the main development 
of this paper. A generalized Detweiler- Whiting axiom like the one derived in section 
[4] will trivially exist for an effective momentum 

V' £ = V^ + £' e , (A A) 

where 



JE L Jo 



J a dS a . (A.5) 



These definitions are easily applied to a family charge distributions defined by (|4"4")) 
and (|45|) . Using the results of section [3] together with multipole expansions for the 
force and torque derived in [171 1 1 8] , it may shown that 



K = qF^ri" - \Q hc ^aFt? + 0(A 3 ), (A.6) 

N' ab = 2Q c la F$ + 0(A 3 ). (A.7) 

These quantities are defined by analogy to ([25")) , where the momenta p' a and S' ab are 
associated with V'^ via a relation like 1(5)). The homogeneous field F^ b m has been 
defined by (j4"Tj) . Note that terms involving the quadrupole and higher moments of J a 
and T ab are negligible in the approximations used here. It may also be shown that 
the 4-velocity 7 never appears anywhere except in the Lorentz force term already 
written down [TS1 [Mj ■ This is essential to the use of a center-of-mass condition as part 
of a well-posed initial value problem. 

For reference, the dipole moment of the current distribution Q ab = Q[ w(s) is 
given by @3] 



J V [a a (e b]c Y J b ' ~ V b] V a ,aJ a 't b '^ dSv, (A.. 



where 



Mx, 7 ) = 0(«&) = / V a V Q 'CrVf,V a 'crdK. (A.9) 

Jo 

As in (|A.1[) , this integral is to be carried out along the geodesic connecting 7 and 
x. Note that 9 ab — * g ab in flat spacetime. This remains approximately true for 
all spacetimes in regions regions where x is sufficiently close to 7. Given some 
preferred timelike vector field, Q ab may be decomposed into electric and magnetic 
dipole moments in the same way that F ab can be written in terms of electric and 
magnetic fields. 
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